Autozygome-guided exome-first study in a consanguineous cohort with early-onset retinal disease uncovers an isolated RIMS2 phenotype and a retina-enriched RIMS2 isoform.

Leber congenital amaurosis (LCA) and early-onset retinal degeneration (EORD) are inherited retinal diseases (IRD) characterized by early-onset vision impairment. Herein, we studied 15 Saudi families by whole exome sequencing (WES) and run-of-homozygosity (ROH) detection via AutoMap in 12/15 consanguineous families. This revealed (likely) pathogenic variants in 11/15 families (73%). A potential founder variant was found in RPGRIP1. Homozygous pathogenic variants were identified in known IRD genes (ATF6, CRB1, CABP4, RDH12, RIMS2, RPGRIP1, SPATA7). We established genotype-driven clinical reclassifications for ATF6, CABP4, and RIMS2. Specifically, we observed isolated IRD in the individual with the novel RIMS2 variant, and we found a retina-enriched RIMS2 isoform conserved but not annotated in mouse. The latter illustrates potential different phenotypic consequences of pathogenic variants depending on the particular tissue/cell-type specific isoforms they affect. Lastly, a compound heterozygous genotype in GUCY2D in one non-consanguineous family was demonstrated, and homozygous variants in novel candidate genes ATG2B and RUFY3 were found in the two remaining consanguineous families. Reporting these genes will allow to validate them in other IRD cohorts. Finally, the missing heritability of the two unsolved IRD cases may be attributed to variants in non-coding regions or structural variants that remained undetected, warranting future WGS studies.

Nigral transcriptomic profiles in Engrailed-1 hemizygous mouse models of Parkinson's disease reveal upregulation of oxidative phosphorylation-related genes associated with delayed dopaminergic neurodegeneration.

Introduction: Parkinson's disease (PD) is the second most common neurodegenerative disorder, increasing both in terms of prevalence and incidence. To date, only symptomatic treatment is available, highlighting the need to increase knowledge on disease etiology in order to develop new therapeutic strategies. Hemizygosity for the gene Engrailed-1 (En1), encoding a conserved transcription factor essential for the programming, survival, and maintenance of midbrain dopaminergic neurons, leads to progressive nigrostriatal degeneration, motor impairment and depressive-like behavior in SwissOF1 (OF1-En1+/-). The neurodegenerative phenotype is, however, absent in C57Bl/6j (C57-En1+/-) mice. En1+/- mice are thus highly relevant tools to identify genetic factors underlying PD susceptibility. Methods: Transcriptome profiles were defined by RNAseq in microdissected substantia nigra from 1-week old OF1, OF1- En1+/-, C57 and C57- En1+/- male mice. Differentially expressed genes (DEGs) were analyzed for functional enrichment. Neurodegeneration was assessed in 4- and 16-week old mice by histology. Results: Nigrostriatal neurodegeneration was manifested in OF1- En1+/- mice by increased dopaminergic striatal axonal swellings from 4 to 16 weeks and decreased number of dopaminergic neurons in the SNpc at 16 weeks compared to OF1. In contrast, C57- En1+/- mice had no significant increase in axonal swellings or cell loss in SNpc at 16 weeks. Transcriptomic analyses identified 198 DEGs between OF1- En1+/- and OF1 mice but only 52 DEGs between C57- En1+/- and C57 mice. Enrichment analysis of DEGs revealed that the neuroprotective phenotype of C57- En1+/- mice was associated with a higher expression of oxidative phosphorylation-related genes compared to both C57 and OF1- En1+/- mice. Discussion: Our results suggest that increased expression of genes encoding mitochondrial proteins before the onset of neurodegeneration is associated with increased resistance to PD-like nigrostriatal neurodegeneration. This highlights the importance of genetic background in PD models, how different strains can be used to model clinical and sub-clinical pathologies and provides insights to gene expression mechanisms associated with PD susceptibility and progression.

Multi-omics analysis in human retina uncovers ultraconserved cis-regulatory elements at rare eye disease loci.

Cross-species genome comparisons have revealed a substantial number of ultraconserved non-coding elements (UCNEs). Several of these elements have proved to be essential tissue- and cell type-specific cis-regulators of developmental gene expression. Here, we characterize a set of UCNEs as candidate CREs (cCREs) during retinal development and evaluate the contribution of their genomic variation to rare eye diseases, for which pathogenic non-coding variants are emerging. Integration of bulk and single-cell retinal multi-omics data reveals 594 genes under potential cis-regulatory control of UCNEs, of which 45 are implicated in rare eye disease. Mining of candidate cis-regulatory UCNEs in WGS data derived from the rare eye disease cohort of Genomics England reveals 178 ultrarare variants within 84 UCNEs associated with 29 disease genes. Overall, we provide a comprehensive annotation of ultraconserved non-coding regions acting as cCREs during retinal development which can be targets of non-coding variation underlying rare eye diseases.

Combining a prioritization strategy and functional studies nominates 5'UTR variants underlying inherited retinal disease.

BACKGROUND: 5' untranslated regions (5'UTRs) are essential modulators of protein translation. Predicting the impact of 5'UTR variants is challenging and rarely performed in routine diagnostics. Here, we present a combined approach of a comprehensive prioritization strategy and functional assays to evaluate 5'UTR variation in two large cohorts of patients with inherited retinal diseases (IRDs). METHODS: We performed an isoform-level re-analysis of retinal RNA-seq data to identify the protein-coding transcripts of 378 IRD genes with highest expression in retina. We evaluated the coverage of their 5'UTRs by different whole exome sequencing (WES) kits. The selected 5'UTRs were analyzed in whole genome sequencing (WGS) and WES data from IRD sub-cohorts from the 100,000 Genomes Project (n = 2397 WGS) and an in-house database (n = 1682 WES), respectively. Identified variants were annotated for 5'UTR-relevant features and classified into seven categories based on their predicted functional consequence. We developed a variant prioritization strategy by integrating population frequency, specific criteria for each category, and family and phenotypic data. A selection of candidate variants underwent functional validation using diverse approaches. RESULTS: Isoform-level re-quantification of retinal gene expression revealed 76 IRD genes with a non-canonical retina-enriched isoform, of which 20 display a fully distinct 5'UTR compared to that of their canonical isoform. Depending on the probe design, 3-20% of IRD genes have 5'UTRs fully captured by WES. After analyzing these regions in both cohorts, we prioritized 11 (likely) pathogenic variants in 10 genes (ARL3, MERTK, NDP, NMNAT1, NPHP4, PAX6, PRPF31, PRPF4, RDH12, RD3), of which 7 were novel. Functional analyses further supported the pathogenicity of three variants. Mis-splicing was demonstrated for the PRPF31:c.-9+1G>T variant. The MERTK:c.-125G>A variant, overlapping a transcriptional start site, was shown to significantly reduce both luciferase mRNA levels and activity. The RDH12:c.-123C>T variant was found in cis with the hypomorphic RDH12:c.701G>A (p.Arg234His) variant in 11 patients. This 5'UTR variant, predicted to introduce an upstream open reading frame, was shown to result in reduced RDH12 protein but unaltered mRNA levels. CONCLUSIONS: This study demonstrates the importance of 5'UTR variants implicated in IRDs and provides a systematic approach for 5'UTR annotation and validation that is applicable to other inherited diseases.

Loss of ctnnd2b affects neuronal differentiation and behavior in zebrafish.

Delta-catenin (CTNND2) is an adhesive junction associated protein belonging to the family of p120 catenins. The human gene is located on the short arm of chromosome 5, the region deleted in Cri-du-chat syndrome (OMIM #123450). Heterozygous loss of CTNND2 has been linked to a wide spectrum of neurodevelopmental disorders such as autism, schizophrenia, and intellectual disability. Here we studied how heterozygous loss of ctnnd2b affects zebrafish embryonic development, and larvae and adult behavior. First, we observed a disorganization of neuronal subtypes in the developing forebrain, namely the presence of ectopic isl1-expressing cells and a local reduction of GABA-positive neurons in the optic recess region. Next, using time-lapse analysis, we found that the disorganized distribution of is1l-expressing forebrain neurons resulted from an increased specification of Isl1:GFP neurons. Finally, we studied the swimming patterns of both larval and adult heterozygous zebrafish and observed an increased activity compared to wildtype animals. Overall, this data suggests a role for ctnnd2b in the differentiation cascade of neuronal subtypes in specific regions of the vertebrate brain, with repercussions in the animal's behavior.

The MHC class II transactivator modulates seeded alpha-synuclein pathology and dopaminergic neurodegeneration in an in vivo rat model of Parkinson's disease.

BACKGROUND: Abnormal folding, aggregation and spreading of alpha-synuclein (αsyn) is a mechanistic hypothesis for the progressive neuropathology in Parkinson's disease (PD). Spread of αsyn between cells is supported by clinical, neuropathological and experimental evidence. It has been proposed that a pro-inflammatory micro-environment in response to αsyn can promote its aggregation. We have previously shown that allelic differences in the major histocompatibility complex class two transactivator (Mhc2ta) gene, located in the VRA4 locus, alter MHCII expression levels, microglial activation and antigen presentation capacity in rats upon human αsyn over-expression. In addition, Mhc2ta regulated dopaminergic neurodegeneration and the extent of motor impairment. The purpose of this study was to determine whether Mhc2ta regulates αsyn aggregation, propagation and dopaminergic pathology in an αsyn pre-formed fibril (PFF)-seeded in vivo model of PD. METHODS: The DA and DA.VRA4 congenic rat strains share background genome but display differential microglial antigen presenting capacity due to different Mhc2ta alleles in the VRA4 locus. PFFs of human αsyn or BSA solution were injected unilaterally to the striatum of DA and DA.VRA4 rats two weeks after ipsilateral administration of recombinant adeno-associated virus (rAAV) vectors carrying human αsyn or GFP to the substantia nigra pars compacta. Behavioural assessment was performed at 2, 5 and 8 weeks while histological evaluation of αsyn pathology, inflammation and neurodegeneration as well as determination of serum cytokine profiles were performed at 8 weeks. RESULTS: rAAV-mediated expression of human αsyn in nigral dopaminergic neurons combined with striatal PFF administration induced enhanced αsyn pathology in DA.VRA4 compared to DA rats. Mhc2ta thus significantly regulated the seeding, propagation and toxicity of αsyn in vivo. This was reflected in terms of wider extent and anatomical distribution of αsyn inclusions, ranging from striatum to the forebrain, midbrain, hindbrain and cerebellum in DA.VRA4. Compared to DA rats, DA.VRA4 also displayed enhanced motor impairment and dopaminergic neurodegeneration as well as higher levels of the proinflammatory cytokines IL-2 and TNFα in serum. CONCLUSIONS: We conclude that the key regulator of MHCII expression, Mhc2ta, modulates neuroinflammation, αsyn-seeded Lewy-like pathology, dopaminergic neurodegeneration and motor impairment. This makes Mhc2ta and microglial antigen presentation promising therapeutic targets for reducing the progressive neuropathology and clinical manifestations in PD.

Loss of Function of RIMS2 Causes a Syndromic Congenital Cone-Rod Synaptic Disease with Neurodevelopmental and Pancreatic Involvement.

Congenital cone-rod synaptic disorder (CRSD), also known as incomplete congenital stationary night blindness (iCSNB), is a non-progressive inherited retinal disease (IRD) characterized by night blindness, photophobia, and nystagmus, and distinctive electroretinographic features. Here, we report bi-allelic RIMS2 variants in seven CRSD-affected individuals from four unrelated families. Apart from CRSD, neurodevelopmental disease was observed in all affected individuals, and abnormal glucose homeostasis was observed in the eldest affected individual. RIMS2 regulates synaptic membrane exocytosis. Data mining of human adult bulk and single-cell retinal transcriptional datasets revealed predominant expression in rod photoreceptors, and immunostaining demonstrated RIMS2 localization in the human retinal outer plexiform layer, Purkinje cells, and pancreatic islets. Additionally, nonsense variants were shown to result in truncated RIMS2 and decreased insulin secretion in mammalian cells. The identification of a syndromic stationary congenital IRD has a major impact on the differential diagnosis of syndromic congenital IRD, which has previously been exclusively linked with degenerative IRD.